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Abstract 

The critical temperature and the nature of the QCD finite temperature phase transition are 
determined for Nj = 2 dynamical flavors of nonperturbatively improved Wilson fermions. The 
calculations are performed on large lattices with temporal extents Nt = 12, 10 and 8, and lattice 
spacings down to a = 0.075 fm. We find the deconfinement and chiral phase transitions to take 
place at the same temperature. Our results are in broad agreement with a second order phase 
transition in the chiral limit. The critical temperature at the physical quark mass is found to be 
Tc = 174(3) (6) MeV. 

PACS numbers: 11.15.Ha, 12.38.Aw, 12.38. Gc 
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I. INTRODUCTION 



One of the basic questions in finite temperature QCD is: Wfiat is tfie nature of the finite 
temperature phase transition, and at which temperature does it happen? In spite of enor- 
mous computational efforts jl|, the answer to this question has remained controversiaL 
The Wuppertal group 2! finds the deconfining transition, which we identify with the peak of 
the Polyakov-loop susceptibihty, and the chiral transition at widely separated temperatures. 
In contrast, the Brookhaven/Bielefeld collaboration claims both temperatures to coincide. 
Furthermore, the Brookhaven/Bielefeld collaboration js] quotes a transition temperature 
of Tc = 196(3) MeV, while the transition temperatures found by the Wuppertal group [4| 
are = 170(7) MeV for the deconfining transition and = 146(5) MeV for the chiral 
transition. Both groups use rooted staggered fermions, but with different levels of improve- 
ment. The Brookhaven/Bielefeld collaboration uses asqtad as well as p4 fermions, while the 
Wuppertal group employs twice iterated stout smeared links in the fermion matrix. It has 
been argued that the discrepancy is largely due to the rather coarse lattices used by the 
Brookhaven/Bielefeld collaboration. Indeed, inititial calculations of this group were limited 
to lattices of temporal extent A^^ < 6, corresponding to lattice spacings of a > 0.2 fm, while 
the Wuppertal group performed simulations on lattices of extent Nt = 10, 8 and 6, and at- 
tempted a continuum extrapolation. More recently, the Brookhaven/Bielefeld collaboration 
has extended their calculations to lattices of temporal extent A^t = 8 and found that with 
decreasing lattice spacing Tc shifted by 5 — 7 MeV towards smaller values. 

The connection between deconfining and chiral transition has been the subject of sev- 
eral phenomenological considerations. Naively, one would expect the temperature of the 
deconfinement transition to lie below that of the chiral transition, if different at all. This 
turns out to be the case, for example, in the Polyakov-loop extended Nambu-Jona-Lasinio 
model j^. More likely is that both transitions occur at the same temperature, as Polyakov 
loop and chiral condensate mix at finite dynamical quark masses. The consequence would 
be a simultaneous enhancement of both the chiral and Polyakov-loop susceptibilities along 



the transition line 







io|. 



To clarify the issue, an independent investigation of the nature of the finite temperature 
phase transition is needed. In this work we shall perform simulations with Nf = 2 dynamical 
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flavors of nonperturbatively 0(a) improved Wilson fermions, so-called clover fermions, and 
plaquette gauge action on lattices of temporal extent Nt = 12, 10 and 8. This action has been 
successfully employed by the ALPHA, UKQCD, QCDSF and CERN/Rome collaborations 
in calculations at zero temperature, on whose results we can draw. Besides that, clover 
fermions have an exact flavor symmetry, which we consider a big advantage over staggered 
fermions, as the nature of the finite temperature phase transition largely depends on the 
flavor degrees of freedom. A disadvantage though is the lack of chiral symmetry. Preliminary 



results of our work have been reported in 111 ]. 



The paper is organized as follows. In sec. II we present the action, and in sec. Ill we 
introduce the order parameters, which are the main focus of the analysis. The lattice data are 
presented in sec. IV and in the Appendix. In sec. V we recapitulate the equation of state, that 
is expected to describe the thermodynamic properties of QCD with two dynamical flavors 
near the critical point. Relations between second derivatives of the partition function are 
called Maxwell relations. One such relation connects the chiral susceptibility to the derivative 
of the plaquette with respect to mass, and will be derived in sec. VI. In sec. VII we compute 
the transition temperature from the Polyakov-loop susceptibility, the chiral susceptibility 
and the correlator of Polyakov loop and chiral condensate, and fit its dependence on the 
quark mass to the prediction of the equation of state. Finally, in sec. VIII we conclude. 

II. LATTICE SIMULATION 

We consider mass-degenerate sea quarks. The fermionic action for each of the two flavors 
reads 

SF = a^Yl Yl ^^^^ ^('^^^ [^/^ ~ ^] + 

X fl 

- ^'^^^ ^l^^ ~ [^/^ + ^] '^^^ ~ (1) 

- csw ^ ^ ^(a;) (^tiu P,,u{x) ipix) + {112 + ipix) ipix] 
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/5 


csw 


V = Nt 




ro/a 


5.20 


2.0171 


16^ 8 


0.136050(17) 


5.454(58) 


5.20 


2.0171 


24^ 10 


0.136050(17) 


5.454(58) 


5.25 


1.9603 


16^8 


0.136273(7) 


5.880(26) 


5.25 


1.9603 


24^8 


0.136273(7) 


5.880(26) 


5.25 


1.9603 


32^ 12 


0.136273(7) 


5.880(26) 


5.29 


1.9192 


24^ 12 


0.136440(4) 


6.201(25) 



TABLE I: Parameters of the simulation. 



where P^u is the clover-leaf form of the lattice field strength tensor, 



+Ul{x — aO) Uj^{x — afi — aO) Uy{x — afi — aO) U^{x — afi) 
—U^{x) Uj^{x — afi + au) Ul{x — afi) U^{x — afi) 
—Ul{x — au) U^{x — au) U^{x + afi — aO) Uj^{x) 



(2) 



and 



arrir 



(3) 



1 1 1 
, am = 

Kc being the critical value of the hopping parameter. The improvement coefficient csw was 



determined nonperturbatively 



12|. 



We use the highly optimized HMC algorithm of the QCDSF collaboration 13| for updat- 
ing the gauge field. Considerable speedups were obtained by applying mass preconditoning 
a la Hasenbusch 1^ and putting the pseudofermion action on multiple time scales 15 |. 
The algorithm runs about three times faster than an equally optimized RHMC algorithm 



for A'^^ = 2 + 1 dynamical fiavors [16], which was one of the reasons for concentrating on 
Nf = 2 flavors first. The couplings, lattice volumes and lattice spacings covered by our 
simulations are listed in Table [H The scale parameters ro/a have been taken from the zero 
temperature runs of the QCDSF collaboration at the corresponding couplings. They refer 
to the chiral limit k, = Kc- We also list the critical hopping parameters Kc, which we adopted 
from QCDSF as well. (For recent relevant work see 17|.) On each lattice we have performed 



simulations at up to 13 different k values placed around the transition line. The exact values 
are given in Tables [TTJHY] in the Appendix. The temperature of our lattices is given by 

T=^- (4) 
Nta ^ ^ 

The gauge field configurations were generated on the BlueGene/L at KEK, the RSCC 
cluster at RIKEN, the MVS-50k at the Joint Computer Center (Moscow), on the SKIF- 
Chebyshev at Moscow State University, as well as on the Altix at HLRN. 

III. ORDER PARAMETER 

Two-flavor QCD is expected to undergo a second order transition at finite temperature 
in the chiral limit and at very small quark masses. In the chiral limit the order parameter 
is the subtracted scalar density cr_ = a — ctq, defined through 

V 



(T = — ^i){x)i){x) = Z'^aR + (To , (5) 



where ctr is the renormalized density, and o"o is an additive renormalization constant, which 
arises from mixing of a with the vacuum due to lack of chiral symmetry. For heavy quark 
masses close to the quenched limit, the theory is known to undergo a first order phase 
transition. In that limit the order parameter is the Polyakov loop 

1 1 

^ = ]^ E ^^^^(^) ' ^(^) = 3 n ■ 

* X X4 = l 

In the intermediate mass region we expect the transition to be a crossover. This scenario is 
sketched in Fig. [1] An 'order parameter', which interpolates between the two limits, is 

1 am ^ 

^= . . + (7) 



^yl + {amy [am^ 



The temperature of the chiral transition is, for general m, identified with the peak position 
of the chiral susceptibility 

X. ^ W\ = {a') - {a)' , (8) 
while the peak of the Polyakov-loop susceptibility 

XL^NUL%AL')c=m-{Lf) (9) 



00 



FIG. 1: The phase diagram of two-flavor QCD. The solid lines indicate a second (bottom) and first 
order transition (top), respectively, while the dashed line denotes the crossover region. 

defines the temperature of the deconfining transition. Both, (cr^)c and {L'^)c can be combined 
into one susceptibility, following ([7]), 

1 , o, 2am , {amy 



[am] 



1 + (am)2 



1 + {am 



where 



(La), = (La) - (L) (a) 



(10) 



(11) 



which interpolates between zero and infinite quark mass. In case the crossover temperature 
is unique, all three correlators, (L^),, (La), and (a^),, are expected to peak at the same 
transition temperature. If, on the other hand, the deconfining and chiral transitions take 
place at far different temperatures, we would not expect to find a distinct peak in the 
correlator (La),. In the connected correlators (a^), and (La), the additive renormalization 
constant ag drops out, which allowed us to replace a_ by a. 

IV. DATA 



To compute the chiral susceptibility Xa ^iid the correlator of L and a, {La)c, all we need 
to know is the average plaquette 

P = lTr ([/□), (12) 
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FIG. 2: The average plaquette (O) on the 24^ 10 (top) and 16^ 8 lattice (bottom) at ^ = 5.20, 
together with the average plaquette at zero temperature from the 24^48 lattice (□). The dashed 
line indicates the position of Kc- 
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FIG. 3: The average Polyakov loop on the 24^ 10 (top) and 32^ 12 lattice (bottom) at /3 = 5.20 
and 5.25, respectively. The solid line in the bottom figure denotes the integral of the Gaussian in 
Fig. El (bottom). 
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FIG. 4: The Polyakov-loop susceptibility on the 16^8 (top), 24^10 (middle) and 32^12 lattice 
(bottom) at (3 = 5.20, 5.20 and 5.25, respectively, together with a Gaussian adaptation. 
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and the average Polyakov loop (L), as we shall see. In Tables UTTI - IVl in the Appendix we 
present our results for P and (L), as well as for the Polyakov-loop susceptibility xl- The 
statistical errors are computed by blocked jackknife. In Fig. |2]we plot the average plaquette 
P on the 24^ 10 and 16^8 lattice at /3 = 5.20. For comparison, we also show the respective 
numbers at zero temperature, taken from the QCDSF collaboration. Both sets of numbers 
agree with each other at smaller k values (corresponding to the confined phase - see Fig. [5]), 
as expected, while at larger k values the finite temperature plaquette exceeds the zero 
temperature one, moving closer to its perturbative value. In Fig. |3] we plot the average 
Polyakov loop (L) on the 24^ 10 and 32'^ 12 lattice at /3 = 5.20 and 5.25, respectively. In 
both cases (L) shows a sharp increase in a narrow interval of k. Finally, in Fig. H] we show 
the Polyakov-loop susceptibility for three different volumes and A^t's. While we observe a 
distinct peak on the 16'^ 8 lattice, it appears that the smaller the quark mass is, the more is 



the peak washed out. This phenomena has also been observed by the CP-PACS group |18 |. 
(Unfortunately, the Wuppertal group does not show the Polyakov-loop susceptibility, so we 
cannot compare.) 

V. EQUATION OF STATE 

We shall assume for the moment that the finite temperature transition of two-flavor QCD 
is indeed of second order in the chiral limit, as indicated in Fig. [H Later on we shall see that 
our data support this view, while a first order transition appears to be unlikely. Throughout 
this section we shall use a and rh as shorthand for the subtracted chiral condensate 
and am, respectively. 

In the vicinity of the phase transition the chiral condensate, the dynamical quark mass 
and the temperature are related by the equation of state 

m = a^f{t/m7^), (13) 

where S and (3 are the critical exponents characterizing the transition, and 

Um = 0) ^'^^ 
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is the reduced temperature, with Tc(m = 0) being the critical temperature in the chiral 



0.537(7) 
0.206(1) 



(15) 



hmit. It is expected that the two-flavor theory is in the same universality class as the 
three-dimensional 0(4) Heisenberg model [l^, with the external magnetic field and the 
magnetization being identified with the bare quark mass m and the chiral condensate a, 
respectively. The critical exponents of this model were found to be 201] 

1 

1 

6 

We may expect 0(4) scaling to be realized in two-fiavor lattice QCD only if the fermion 
action obeys fiavor and chiral symmetry. Neither of the common actions satisfy both. The 
problem of 0(4) scaling has been investigated by several authors (for a review see 21|) with 
somewhat contradictory results. Simulations of two-fiavor rooted staggered fermions 22] 
show deviations from 0(4) scaling, albeit at small Nt. Since staggered fermions show quite 
large taste breaking effects at strong coupling, one mip:ht expect to find 0(2) critical expo- 
nents instead. This conjecture receives support from The situation improves with the 
adoption of p4 fermions. In this case evidence for 0(N) scaling was found just recently 24 1. 
For both Wilson 25|] and clover fermions 26|] , on the other hand, it has been shown that the 
subtracted condensate scales with 0(4) critical exponents (for Tc(m) > Tc(0) though, see 
below). In 26|] good agreement was found even for pion masses up to 1 GeV. At our lattice 
spacings, a ^ 0.08 fm, we find chiral symmetry to be largely restored, which manifests itself, 



for example, in the equa! 

Zp 



J •! • ry Singlet 

densities, ^ 



ity of the renormalization constants of the scalar and pseudoscalar 



271 ]. So we may hope for 0(4) scaling. 



Let Tc(m) denote the pseudocritical temperature at finite m, which we define to be the 
temperature corresponding to the position of the peak of the chiral susceptibility 

da 



dm 



From the scaling relation f ll3p we then derive 

Tc{m) — Tc{m = 0) oc . 



(16) 



(17) 



Assuming 



m_ oc m 



(18) 
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we thus expect to find 

Te(m) - T,(m = 0) oc m^f^^^^ (19) 

for a second order transition at m = 0. A first order transition, on the other hand, would 
give 



Tc{m) — Tc{m = 0) oc . 



Ahernatively, at fixed Nf, the reduced temperature t can be identified with 25 1 

t = /3-f3 



(20) 



(21) 



where is the value of (3, at which the critical temperature Tc{m) line hits the critical 
hopping parameter line as we lower (3. In Fig. [5] we show the line of Tc{m) in the k — (3 
plane for Nt = 8, 10 and 12. 

Expanding the scaling function / in ( ITSl) to lowest nontrivial order in t, i.e. neglecting 
terms of O(t^), we obtain 

rh = Aa^ + Bta^~^ (22) 

with A,B>0. From (!22|) we can immediately read off the phase diagram in the k — 13 plane. 
This is shown in Fig. [6] for Nt = 12. Similar diagrams hold for different values of Nt. For 
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FIG. 5: The lines of critical temperature Tc{m) for Nt = 8, 10 and 12, together with the critical 
hopping parameter Kc line. The open circles refer to Tc{m) obtained from the chiral susceptibility. 
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FIG. 6: The phase diagram of clover fermions for Nt = 12. 

t > the chiral condensate vanishes in the chiral hmit. For t < 0, on the other hand, chiral 
symmetry is spontaneously broken for all k values and 

B 



a/3 



A 



t > 



(23) 



in the chiral limit. The critical temperature Tc(m = 0) is obtained from Nt and the corre- 
sponding lattice spacing at t = 0. 

VI. CHIRAL SUSCEPTIBILITY AND MAXWELL RELATION 



The chiral condensate is related to the average plaquette P by means of a Maxwell 



relation 



28 



The plaquette can be computed with high precision and so lends itself to an 
accurate determination of the critical temperature of the chiral transition. 

Both the chiral condensate and the plaquette can be found from the partial derivatives 



of the partition function Z: 



V df3 

V dm 



-6P + 2 



drhc 



a -2 



dcsw_ 
dp 



2a. 



(24) 
(25) 



where, temporarily, we have set a = (a), and 



'^=W' <^= 2 ^S1^(^)^'^'^^^''(^)^(^)• 



(26) 
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The second derivative d"^ \n. Z / d (i dm can be expressed in two different orders, 



Vd(5dm ^ dp 



-6 



dP 



dm 



drhf. da 
+ 2 

^ d (5 dm 



d csw 9 6 



d (5 dm 



(27) 



Both 0" and 5 are chiral order parameters. In the following we shall neglect the contribution 
from the clover term 5, mainly for the sake of clarity. We may do so, because it is suppressed 



by two orders of the lattice spacing with respect to the chiral condensate 

5 



29|-|3l|: 



This leaves us with the relation 



dP 



dm 



1 drhc da 



/3 



2> d j3 dm 



^ = 0{a') . 
a 



1 da 
3 



1 dm 



3 dp 



da 



dm 



(28) 



(29) 



Above identity is called the Maxwell relation. It holds for any lattice size and for all values 
of /3 and m. The right-hand form of the identity is the most useful, as it gives us the chiral 
susceptibility: 



da 



dm 



drhc 



+ 



dm 


J 


dP 




^ d 1 


.) 


dP 


df3 




dm 








dm 



(30) 



In both correlators, Xo- a-^d {d P/dm)\^, the additive renormalization ao drops out, so that 
we may identify a with (a^) again. 

Before we can identify any peak of dP/drh with a peak of Xo-, we need to make sure that 
the factor in front of dP/drh, 

-1 



/i 



drhc dm 
+ 



dp dp 

is 'well behaved'. The QCDSF collaboration finds for the first term 



(31) 



drhc 
~dj 



-0.1536 P = 5.20 
-0.1177 /3 = 5.25 . 
-0.0919 P = 5.29 



(32) 



The second term can be derived from the equation of state fl22!) . which gives 



dm 
~dp 



B 



a 



(33) 
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FIG. 7: The chiral susceptibility on the 16^8 (top), 24? 10 (middle) and 32^ 12 lattice (bottom) at 
j3 = 5.20, 5.20 and 5.25, respectively, together with a Gaussian. 
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FIG. 8: The correlator {La)c on the 16^8 (top), 24^10 (middle) and 32^12 lattice (bottom) at 
j3 = 5.20, 5.20 and 5.25, respectively, together with a Gaussian. 
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It vanishes proportional to m in the chiral hmit, so that for small masses fi~^ is dominated 
by the first term fl32|) . 

It does not matter whether we use the renormalized or unrenormalized subtracted scalar 
density to determine the transition temperature, because Zs is an extremely smooth function 
of rh 32|, which varies by less than 1% over the transition region. 



The partial differential equation (129|1 can be solved by 



where we have made use of the fact that a = in the chiral limit. For the ansatz fl33|) we 
find 

B A 1 4.1 1 drhf. , . 

p-p("^-')=w^)' "^-^iff'- '''' 

Knowing the critical exponents and 5, we may compute the chiral condensate from the 
average plaquette, and vice versa. 

The chiral transition has merely a small effect on the average plaquette. The reason is 
that only a fraction of a percent of the magnitude of P is of nonperturbative origin 33 |. 
Nevertheless, with accurate data of the plaquette the Maxwell relation proves to be a viable 
tool to locate the position of the chiral phase transition. We use the symmetric difference 
quotient method, 

/'W = ^'"^ (36) 

to compute the derivative of the plaquette, dP/dm\j3. In Fig.[71we plot the chiral condensate 
jj/^Q- as a function of k for three representative lattices, which we chose to be the same as 
those in Fig. |H 



VII. TRANSITION TEMPERATURE(S) 

To complete the calculation of the susceptibility (w^)c introduced in (fTOl) . what remains 
to be computed is the correlator (Lcr)^. This can be obtained from the derivative of the 
average Polyakov loop with respect to mass: 

d{L) 



^ '^'^^ dm 



(37) 
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FIG. 9: The Polyakov-loop susceptibility xl (^7), the chiral susceptibihty Xa (^4) and the corre- 
lator (L(t)c on the 16^ 8 lattice at /3 = 5.20. 

Again, the additive renormalization constant drops out i n fl37|) . and it does not make any 
difference whether we use the renormahzed Polyakov loop 34] or not, because the renormal- 
ization factor is basically a function of (3 only, being an ultraviolet quantity. In Fig. |8] we 
plot {La)c on the lattices of Figs. [Hand [71 where we have expressed the derivative by the 
symmetric difference ( l36l) as before. Not only do we find that the Polyakov loop and the 
chiral condensate are strongly correlated within a narrow temperature range, but also that 
the position of the peaks of all three quantities, xl, Xcr and {La)c, coincide within small 
error bars. This is illustrated once more in Fig. [9] for the 16^ 8 lattice at /3 = 5.20. 

We may express the k values by the corresponding pion masses at zero temperature. The 
latter are known from simulations of the QCDSF collaboration. In Table [IT] we give the 
pseudocritical temperature Tc(m) and the corresponding pseudocritical pion masses, m^", 
obtained from the peak of the Polyakov-loop susceptibility, the chiral susceptibility and the 
correlator ( 137]) of L and a, respectively, and in Fig. [10] we plot the results together with an 
extrapolation to the chiral limit. On all lattices the individual pion masses m^" are found 
to coincide with each other within the error bars. We do not observe any scaling violations. 
Our complementary runs at zero temperature did not reveal any visible scaling violations 
either. The spatial volume varies from Ls{= aNg) = 1.4 fm to Lg = 2.5 fm, while the aspect 
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/3 


V 




XL 


T 

T*!-! tin 
' '''-Tr 

Xa 


It \ 


5.20 


16^8 


0.682(7) 


2.73(6) 


2.78(6) 


2.81(7) 


5.20 


24^ 10 


0.545(6) 


1.59(8) 


1.59(16) 


1.55(14) 


5.25 


24^8 


0.735(3) 


3.18(4) 


3.17(4) 


3.33(7) 


5.25 


323 12 


0.490(2) 


1.00(11) 


1.05(8) 


1.05(7) 


5.29 


24^12 


0.517(2) 


1.49(8) 


1.40(9) 


1.3(1) 



TABLE II: The pseudocritical temperatures and corresponding pion masses obtained from the 
peak of XL, Xa and {La)c on our various lattices. 

ratio Ng/Nt varies between 2 and 3. Remarkably, all our results fall on a single curve. In 
addition, we do not see any finite size effect in direct comparison at A't = 8. 

The temperature Tc(m) shows an almost linear behavior in the pion mass, in accord with 




FIG. 10: The pseudocritical temperature Tc{m) as a function of pion mass, together with a fit to 
the power ml:^'^ , according to the three-dimensional 0(4) modeL 
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the prediction (fT9l) of the 0(4) modeL We thus may fit the data by the ansatz 

T,{m)=C + D{rom^y-^' . (38) 

The result is shown by the sohd curve. Setting the scale by the nucleon mass, the QCDSF 
collaboration finds tq = 0.467(15) fm. Using this value, we obtain at the physical pion mass 

roTc = 0.412(6) , = 174(3)(6) MeV , (39) 

where the first error on is statistical, and the second error reflects the uncertainty in setting 
the scale. This result is in good agreement with the deconfining transition temperature found 
by the Wuppertal group, but lies significantly below the result of the Brookhaven/Bielefeld 
collaboration. 



VIII. CONCLUSIONS AND OUTLOOK 



We have simulated QCD at finite temperature with two dynamical flavors of nonpertur- 
batively improved Wilson fermions on lattices as large as A^^^ = 12 and lattice spacings as 
low as 0.075 fm. The transition temperature has been computed from the Polyakov-loop 
susceptibility, the chiral susceptibility as well as the correlator of Polyakov loop and chiral 
condensate. All three temperatures are found to coincide with each other within the error 

n 

bars. Our results do not support the claim of the Wuppertal group [2|, albeit for A'^y- = 2 + 1 
flavors, that the deconfining and chiral transitions take place at distinctly separated tem- 
peratures. The critical behavior appears to be in accord with the predictions of the 0(4) 
Heisenberg model, at least as far as the quark mass dependence of Tc is concerned, while 
a first order transition |35|] is very unlikely. However, further simulations at smaller quark 
masses are needed in order to confirm this conclusion beyond doubts. 

The Maxwell relation has proven to be a powerful tool in unveiling the phase structure 
of clover fermions. It would be interesting to test (!34l) . and the equation of state (fT3!) itself, 
by direct calculation of the chiral condensate. As we already mentioned, the clover term 5 
is forbidden by chiral symmetry just like a. So there is actually no need to neglect it (as we 
did) if we are just looking for a quantity which peaks at the chiral transition. 

The next step is to extend the simulations to physical quark masses. Such calculations 
require lattices with temporal extent Nt = 14 and larger, given the fact that simulations 
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at couplings below /3 = 5.20 are not feasible. This requires a major increase in computing 
resources. Preliminary investigations at /3 = 5.25 indicate the existence of a transition below 
Kc, resulting in roTc(m) = 0.420(2). This makes us feel confident that the curve in Fig. [TD] 
will not level off towards the staggered result j^. 

Acknowledgment 

We like to thank the computer centers at KEK (under the Large Scale Simulation Pro- 
gram No. 07-14-B), JSSC (Moscow), SC MSU (Moscow), RIKEN and HLRN (Berlin and 
Hannover) for their generous allocation of computer time and technical support. Further- 
more, we like to thank Thomas Streuer and Hinnerk Stiiben for help with the HMC code. 
VGB, SMM and MIP are supported in part through grants RFBR 08-02-00661, RFBR 09- 
02-0033, DFG-RFBR 436 RUS and NSH-679.2008.2, and like to thank their colleagues at 
RFNC-VNIIEF for invaluable help with the computations. 



22 



APPENDIX: TABLES OF RESULTS 



f3 = 5.20 


16^8 


24^ 10 




P 


(L) 


XL 


K 


P 


(L) 


XL 


0.1330 




0.0016(3) 


0.075(3) 


0.1348 


0.467177(33) 


0.0017(2) 


0.066(4) 


0.1335 




0.0026(6) 


0.095(8) 


0.1352 


0.464857(63) 


0.0027(2) 


0.065(4) 


0.1340 




0.0042(6) 


0.099(8) 


0.1353 


0.464367(38) 


0.0033(3) 


0.074(6) 


0.1343 


0.469500(166) 


0.0065(6) 


0.136(12) 


0.1354 


0.463854(38) 


0.0035(2) 


0.080(4) 


0.1344 


0.468949(189) 


0.0096(7) 


0.153(11) 


0.1355 


0.463204(47) 


0.0057(2) 


0.094(5) 


0.1345 


0.468436(116) 


0.0102(7) 


0.167(15) 


0.1356 


0.462574(42) 


0.0063(2) 


0.082(4) 


0.1346 


0.467274(118) 


0.0147(6) 


0.188(13) 


0.1357 


0.462027(35) 


0.0070(2) 


0.082(3) 


0.1347 


0.466334(102) 


0.0172(8) 


0.182(13) 


0.1358 


0.461536(24) 


0.0079(2) 


0.078(3) 


0.1348 


0.465717(101) 


0.0195(6) 


0.160(10) 


0.1359 


0.461167(37) 


0.0081(2) 


0.075(4) 


0.1349 


0.465274(98) 


0.0198(6) 


0.149(13) 


0.1360 


0.460670(36) 


0.0088(2) 


0.077(6) 


0.1355 




0.0293(6) 


0.111(8) 










0.1360 


0.460441(57) 


0.0290(4) 


0.121(7) 











TABLE III: The average plaquette, the average Polyakov loop and the Polyakov-loop susceptibihty 
on the 16'^ 8 and 24^ 10 lattice at /3 = 5.20 against k. The data is based on 0(10, 000) trajectories at 
K values in the immediate vicinity of the peak of the Polyakov-loop susceptibility and on 0(5, 000) 
trajectories towards the edges. The numbers refer to trajectory lengths of unit one. 
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/3 = 5.25 




24^8 


K 


P 


{L) 


XL 


K 


P 


{L) 


XL 


0.1330 


0.46879(10) 


0.0025(7) 


0.07(1) 


0.1340 


0.46326(7) 


0.0090(4) 


0.144(14) 


0.1337 




0.0057(8) 


0.158(20) 


0.1341 


0.46258(6) 


0.0127(4) 


0.160(11) 


0.13375 


0.46454(10) 


0.0065(6) 


0.151(10) 


0.1342 


0.46197(6) 


0.0148(5) 


0.205(22) 


0.1339 


0.46370(11) 


0.0084(7) 


0.173(19) 


0.1343 


0.46151(7) 


0.0164(5) 


0.231(19) 


0.1340 


0.46315(9 ) 


0.0105(7) 


0.170(16) 


0.1344 


0.46080(6) 


0.0188(5) 


0.224(19) 


0.1341 


0.46267(10) 


0.0115(7) 


0.190(18) 


0.1345 


0.46030(4) 


0.0205(3) 


0.199(17) 


0.1342 


0.46188(9) 


0.0145(6) 


0.203(19) 


0.1346 


0.45985(10) 


0.0220(4) 


0.172(16) 


0.13425 


0.46152(8) 


0.0167(7) 


0.229(19) 


32^ 12 


0.1343 


0.46136(8) 


0.0164(5) 


0.207(11) 


K 


P 


{L) 




0.1344 


0.46062(12) 


0.0182(8) 


0.200(14) 


0.1359 


0.45608(2) 


0.00193(6) 


0.0626(13) 


0.1345 


0.46044(26) 


0.0201(7) 


0.221(17) 


0.13595 


0.45591(3) 


0.00197(8) 


0.0613(20) 


0.1350 


0.46039(16) 


0.0267(10) 


0.163(17) 


0.1360 


0.45571(2) 


0.00218(7) 


0.0627(18) 


0.1361 


0.45479(8) 


0.0328(4) 


0.144(8) 


0.13605 


0.45544(2) 


0.00249(7) 


0.0649(17) 










0.1361 


0.45527(2) 


0.00270(6) 


0.0634(23) 










0.13615 


0.45509(3) 


0.00277(8) 


0.0631(20) 



TABLE IV: The average plaquette, the average Polyakov loop and the Polyakov-loop susceptibiUty 
on the 16^ 8, 24^ 8 and 32^ 12 lattice at /3 = 5.25 against k. The data is based on 0(10, 000-15, 000) 
(0(10,000)) trajectories at k values in the immediate vicinity of the peak of the Polyakov-loop 
susceptibility and on 0(2,000 — 4,000) (0(3,000 — 5,000)) trajectories right at the edges on the 
16^ 8 (24^ 8) lattice. On the 32^ 12 lattice we have accumulated 0(3, 000 - 7, 000) trajectories at 
the individual k values. 



24 



(3 = 5.29 243 12 


K 


P 


(L) 


XL 


0.1357 


0.45223(4) 


0.00165(15) 


0.0665(38) 


0.1358 


0.45185(2) 


0.00184(10) 


0.0667(18) 


0.1359 


0.45153(2) 


0.00238(11) 


0.0695(21) 


0.1360 


0.45114(3) 


0.00272(11) 


0.0684(19) 


0.1361 


0.45086(4) 


0.00340(11) 


0.0623(27) 



TABLE V: The average plaquette, the average Polyakov loop and the Polyakov-loop susceptibihty 
on the 24'^ 12 lattice at (3 = 5.29 against k. The data is based on 0(5,000) trajectories at the 
central k values and on 0(2,000 — 3,000) trajectories at the outer k values. 
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